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ABSTRACT 

The extinction law and the variation of colour excess with position, luminosity as well 
as spectral class in young open star clusters NGC 663, NGC 869, NGC 884, NGC 
1502, NGC 1893, NGC 2244, NGC 2264, NGC 6611, Tr 14, Tr 15, Tr 16, CoU 228, Tr 
37 and Be 86 have been studied. The difference in the minimum and maximum values 
of E{B — V) of cluster members has been considered as a measure of the presence of 
non-uniform gas and dust inside the clusters. Its value ranges from 0.22 to 1.03 mag 
in clusters under study, which indicates that non-uniform extinction is present in all 
the clusters. It has been noticed for the first time in NGC 1502 and Tr 37. It is also 
found that the differential colour excess in open clusters, which may be due to the 
presence of gas and dust, decreases systematically with the age of clusters indicating 
that matter is used either in star formation or blown away by hot stars or both. There 
is no uniformity in the variation of E{B — V) with either position or spectral class or 
luminosity. Except in Tr 14, all clusters show a random spatial distribution of E{B — V) 
indicating a random distribution of gas and dust inside the clusters. The E(B — V) 
value correlates with both luminosity and spectral class only in the case of Coll 228, 
Tr 16 and Be 86. The members of these clusters at A > A^, show larger values of colour 
excess ratios than the normal ones. The value of E{U — V) /E{B — V) for most of the 
cluster members is close to the normal interstellar value of 1.73. However, the colour 
excess ratios with E{B — V) at A>Aj are smaller than the normal value for NGC 663, 
NGC 869, NGC 884 and NGC 1502 while they are larger for NGC 6611, CoU 228, Tr 
16 and Tr 14. Thus there is no uniformity in the relationship of extinction properties 
amongst the clusters under study. 
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1 INTRODUCTION 

The investigation of interstellar matter in star-forming envi- 
ronments is a front line research area. One of the most obvi- 
ous properties of interstellar matter that can be determined 
observationally is the extinction law, i.e. the dependence of 
the extinction properties of the matter on wavelength. This 
has been particularly fruitful in recent years with the exten- 
sion of observing capabilities beyond the optical, into the 
ultra-violet (UV) and infrared. 

In the pre-main sequence evolutionary scenario sug- 
gested by Larson (1973) and Hayashi (1970) a star forms 
in a cloud and contracts gravitationally until a sufficiently 
high central temperature, density and pressure are reached 
for nuclear reactions to begin. As this stellar core evolves to 
main sequence (MS) stability, it continues to accrete the sur- 
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rounding circumstellar envelope of gas and dust. Obviously, 
not all the circumstellar material becomes part of the star; 
some fraction is left as a remnant envelope of gas and dust. 
This view is supported by the infrared observations of pre- 
MS objects in NGC 2264 by Warner, Strom & Strom (1979) 
and in the Orion nebula cluster by McNamara (1976). Also, 
not all matter of the molecular cloud is used in forming stars. 
The young (age < lO'^ years) star clusters thus provide an 
ideal opportunity for studying the properties of matter in 
difi'erent star forming regions containing hot O and B stars. 

The non-uniform extinction in open star clusters of in- 
termediate age (< 10* yrs) was studied by Burki (1975) 
using photoelectric UBV photometric values of early type 
cluster members. He found that differential extinction is 
important not only in the case of very young (~ 10^ yrs) 
clusters, but also in some clusters older than 5 x 10^ years. 
Wallenquist (1975) studied non-uniform extinction in open 
clusters using star counts from the Palomar Sky Survey. 
The underlying assumption in this type of investigation is 
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Table 1. General information taken from Mermilliod (1995) about tlie clusters under study. The Galacto- 
centric distance, no. of stars used in this analysis, statistically expected number of field stars in the sample 
and median proper motion membership probability of the sample are denoted by Rq, Ns, Nf and Pm 
respectively. Rq is calculated assuming galacto-centric distance of the sun as 8.5 kpc. In the source column 
1, 2, 3, 4 and 5 denote references Muminov (1983), Cudworth et al. (1993), Marschall, van Altena & Chiu 
(1982), Zhao et al. (1985) and Marschall & van Altena (1987 ) respectively 
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that the observed deficiency of stars is mainly due to the 
presence of absorbing matter between the observer and the 
cluster. Sagar (1987) studied interstellar extinction in 15 
open clusters using UBV photoelectric photometric obser- 
vations of proper motion cluster members, and found that 
ten of them show non-uniform extinction across the cluster 
region. Reddish (1967) has shown that reddening increases 
with stellar luminosity in all clusters and associations with 
ages < 10^ years, but only in some objects with ages be- 
tween 10^ and 2 x 10^ years, and such relations are not 
observed in the clusters older than 2 x 10® years. Bohan- 
nan (1975) re-examined the young clusters data of Reddish 
(1967) and after identifying foreground stars and applying 
correct intrinsic colour indices for bright supergiants found 
no correlation of reddening with luminosity. Sagar (1987) 
observed that in some young clusters the variation of colour 
excess E{B — V) correlated with luminosity, in the sense that 
brighter cluster members were more highly reddened. Most 
of the above analyses were based on UBV data taken primar- 
ily with either photographic plates or single-channel pho- 
tometers. Current optical and infrared imaging technology 
allows us to derive more accurate photometry, as the effects 
of nebular surface brightness which are generally present in 
young clusters can be properly removed. Recently, such ob- 
servations have become available in the literature. We have 
therefore used them to study the non-uniform extinction as 
well as its nature in 14 young open clusters. Section 2 de- 
scribes the selection of the sample. In Section 3 the details 
of the reddening determinations are given, while the results 
derived from the present analysis and their discussions are 
given in the remaining part of the paper. 



2 SELECTION OF CLUSTERS 

This section provides information about the criteria adopted 
for selection of clusters and their members along with infor- 
mation about observational data. 



2.1 Star clusters and observational data 

The observations used in this study are taken from the com- 
pilation database on star clusters by Mermillio d (1995) at 
web site |itt p : / /obswww ■ unige . ch / webda / , We have se- 
lected those 14 young (age < 20 Myr) open clusters which 
have JHK data for at least 10 cluster members. The general 
information about them are listed in Table 1. The distances 
of the sample clusters range from 0.8 to 3.9 kpc. They are 
distributed non-uniformly along the galactic plane with lon- 
gitude ranging from about 15 to 290 deg. The cluster size 
ranges from ~ 1 to 10 pc with an average radius of ~ 4.5 pc, 
while the galacto-centric distance ranges from ~ 6 to 12 kpc. 
Eight clusters of the sample have UBVJHK data while re- 
maining six have UBVRIJHK data. In those cases in which 
the RI photometry was reported in the Kron-Cousin system, 
the colours are converted to Johnson's VRI system using the 
relations given by Bessell (1979). We have used homogeneous 
and accurate data if they are available from more than one 
sources. The observations used in this present analysis are 
mostly based on modern optical and infrared imaging ex- 
cept for NGC 1502, Tr 15 and Tr 37. They have accuracies 
generally better than 0.02 mag in V, (B— V), (V— R) and 
(V— I) and 0.03 mag in (U— B). The accuracy of JHK data 
are generally ~ 0.07 mag. The spectroscopic data wherever 
available are also taken. The accuracy of MK classification 
is generally better than two subclasses in spectral type and 
a class in luminosity. 



2.2 Cluster Membership 

Selection of cluster members is important for the study of 
interstellar extinction across the cluster region. Among the 
most commonly used methods, the one based on kinemat- 
ical data (proper motion and radial velocity) is considered 
the most reliable. As reliable proper motion data are avail- 
able for 8 of the 14 clusters under study, we use them to 
select cluster members. The histograms of the membership 
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Figure 1. Proper motion membership probability distribution of 
stars in the regions of NGC 869, NGC 884, Tr 14, Tr 16 and Tr 
37. The source of data is indicated in Table 1. 



probability (P) of all the investigated stars in clusters NGC 
869, NGC 884, Tr 14, Tr 16 and Tr 37 are shown in Fig. 
1. For NGC 2244, NGC 2264 and NGC 6611 such diagrams 
are available in the sources listed in Table 1. The histograms 
of P distribution indicate that in a cluster, stars are gener- 
ally either in the low P group ( < 20%) or high P group 
(> 50%). It means that the proper motion studies of these 
regions have been able to segregate cluster members from 
field stars. The probability of inclusion of field stars in the 
sample is reduced if stars with P> 50% are considered as 
cluster members and the same has been done here. For the 
five clusters NGC 1502, NGC 1893, CoU 228, Tr 15 and Be 86 
the photometric criteria have been used to segregate cluster 
members from field stars. For NGC 663, photometric crite- 
ria has been applied to the sample of stars having P>1%, 
as the proper motion data alone are unable to separate the 
cluster members from field stars. 



3 DETERMINATION OF COLOUR EXCESSES 

To study the behaviour of the reddening and extinction law 
across the cluster region, we derive E{U — V), E{B — V), 
E{V-R), E{V-I), E{V-J), E{V-H), E{V - K) colour 
excesses of the cluster members. This was accomplished by 
comparing the observed colours of the stars with their intrin- 
sic colours derived from the MKK spectral type-luminosity 
class colour relation given by FitzGerald (1970) for (U— V) 
and (B-V); by Johnson (1966) for (V-R) and (V-I); and 
by Koornneef (1983) for (V-J), (V-H) and (V-K). MK 
classification based on relatively higher dispersion was pre- 



ferred in cases where it is available from more than one 
sources. In order to have a statistically significant number 
of sample stars, photometric spectral types are determined 
using (U— B) and (B— V) data and the Johnson & Morgan 
(1953) Q - method (see Sagar & Joshi 1979 for details) for 
those stars whose positions in the (U— B), (B— V) and V, 
(B— V) diagrams of the clusters (given in the literature) in- 
dicate that they are early type MS stars. In the clusters, the 
number of stars used in the present analysis are given in Ta- 
ble 1. As the median P values of the cluster members range 
from about 70 to 95%, expected number of field stars in the 
sample range from 5 to 30% (see Table 1). The presence of 
some non-members in the sample may therefore not affect 
the conclusions of the present analysis. Wherever possible, 
homogeneous photometric data have been used. All these 
were done to increase the accuracy of colour excess estima- 
tion. E{U - V), E{B - V), E{V - R) and E{V - I) are gen- 
erally uncertain by ~0.08 mag, while E{V - J), E{V ~ H) 
and E{V — K) have typical errors of ~0.15 mag. This is 
mainly because UBVRI magnitudes are more accurate than 
JHK magnitudes. 



4 PRESENCE OF NON-UNIFORM 
EXTINCTION 

To see the extent of non-uniform extinction in different clus- 
ters under study, we plot the histograms of E{B — V) in 
Fig. 2. This indicates that clusters have a different amount 
of non-uniform reddening as there is a wide range in the 
E{B — V) values of cluster members. 

To see the presence of non-uniform extinction in a clus- 
ter, we calculate the value of AE{B — V) = E{B — V)max 
- E{B - V)min, where E{B - V)max and E{B - V)min 
are determined on the basis of, respectively, the five highest 
and five lowest E{B — V) values of the MS cluster mem- 
bers. The values of AE{B — V) obtained in this way along 
with other relevant informations are listed in Table 2. Apart 
from the non-uniform extinction, the other factors for the 
observed dispersion in E{B — V) are stellar evolution, stel- 
lar duplicity, stellar rotation, difference in chemical compo- 
sition, dispersion in ages, dispersion in distances, presence of 
non member stars and inaccuracies in the photometric data 
(cf. Burki 1975, Sagar 1987). However, they (excluding the 
differential extinction) can produce a maximum dispersion 
in AE{B - V) of about 0.0 - 0.11 mag for MS stars. So, a 
value AE{B ~ V)> 0.11 mag for a cluster is considered as 
an indication of the presence of non-uniform gas and dust 
in the direction of the clusters. This conclusion is also sup- 
ported by near-infrared photometric studies of young open 
star clusters (cf. Tapia et al. 1988; Roth 1983; Sagar & Qian 
1989). Since AE{B — V) values for all the clusters under 
study range from 0.22 to 1.03 mag, they indicate the pres- 
ence of non-uniform extinction in the direction of the clus- 
ters. We thus confirm the results obtained earlier by other 
investigators, e.g. Tapia et al. (1991) in NGC 663; Tapia et 
al. (1984) in NGC 869 and NGC 884; Vallenari et al. (1999) 
in NGC 1893 and Be 86; Perez, The & Westerlund (1987) in 
NGC 2244 and NGC 2264; Sagar (1987) and Hillenbrand et 
al. (1994) in NGC 6611 and Tapia et al. (1988) in Tr 14, Tr 
15, Tr 16 and Coll 228. The presence of non-uniform extinc- 
tion has been indicated in NGC 1502 and Tr 37 for the first 
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Table 2. The area investigated around central coordinates of the clusters under study are listed. The values 
of mean E{B - V), E{B - V)min, E{B - V)max and /^E(B - V) are also given. 
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Figure 2. The presence of non-uniform extinction in all clusters 
under study indicated as the dispersions in E{B — V) are signifi- 
cantly larger than that expected from other sources (see text). 



time. In order to understand the possible reasons for non- 
uniform extinction in the young star clusters, we carried out 
the following studies 



4.1 Variation of AE{B ~ V) with age 

Young star clusters (age < 10^ yr) are embedded in dust 
and gas clouds and it is assumed that with time, gas and 
dust in these clouds will either be used up in star formation 
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Figure 3. Variation of AE{B - V) (= E{B - V)max - E{B - 
V)min) with cluster age. The values of AE{B — V) are taken 
from different sources as indicated in the plot. The open circle 
and associated error bar denote the mean and standard deviation 
of the binned data (see text). 



processes or will be dispersed by radiation pressure due to 
massive stars present in these systems or both. Therefore 
in young open clusters, the presence of a variable amount 
of unused gas and dust is expected to cause non-uniform 
extinction. In order to verify this, we plot in Fig. 3 cluster 
age against AE{B — V) which is considered as representa- 
tive of non uniform extinction. The values of AE{B — V) 
for the clusters are taken from the present work (see Ta- 
ble 2); Sagar (1985); Sagar, Munari & de Boer (2001) and 
Pandey, Mahra &i Sagar (1990) while the corresponding ages 
are taken from the respective sources if available. Other- 
wise, they are adopted from Mermilliod (1995). In order to 
see the general variation of mean AE{B — V) with age, we 
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binned the data in various age groups and plot their mean 
AE{B — V) and standard deviation in Fig. 3. This shows 
that there is a variation of AE{B — V) with age of the clus- 
ter. In clusters older than 10* years, there seems to be no 
signature of the presence of gas and dust inside the clusters, 
as the scatter is minimal. A large scatter in AE{B — V) is 
found in clusters younger than 10 Myr indicating the pres- 
ence of varying amounts of gas and dust inside them. The 
observed relation indicates that with time, primordial gas 
and dust are either being used in star formation or being 
blown away by the hot stars present in the clusters, or both, 
as expected. 

4.2 Spatial variation of colour excess E{B — V) 

To study the spatial variation of colour excess across the 
cluster region, we divide the cluster field into equal areas of 
small boxes of size 5'x5' for the clusters NGC 663, NGC 869, 
NGC 884, NGC 1893, Tr 14, Tr 15 and Be 86 of the sample. 
For the clusters NGC 2244, NGC 2264, and NGC 6611 such 
a study has already been carried out by Ogura & Ishida 
(1981), Sagar & Joshi (1983) and Sagar & Joshi (1979) re- 
spectively. The spatial variation of E{B — V) could not be 
studied in Coll 228 as the cluster members occupy a small 
area on the sky while in NGC 1502, Tr 15 and Tr 37, statis- 
tically insignificant number (see Table 1) of cluster members 
denied such a study. Except for Tr 14, there is no system- 
atic correlation of colour excess with position implying that 
gas and dust responsible for variable extinction may be dis- 
tributed randomly within the cluster. Similar behaviour has 
also been found by Sagar (1987) in young clusters NGC 654, 
NGC 2264, NGC 6823, NGC 6913, IC 1805, NGC 6530 and 
NGC 6611. In Tr 14, E{B - V) appears to vary systemat- 
ically with position (see Table 3). The reddening increases 
from east to west in the southern part of the cluster. The 
clusters NGC 663, NGC 869, NGC 884, NGC 1893, Tr 15 
and Be 86 do not show any positional variation of E(B-V). 

4.3 Variation of colour excess E{B — V) with 
luminosity 

In order to study the variation of colour excess with lu- 
minosity, we have plotted colour excess E{B — V) against 
absolute magnitude My for the clusters under study in Fig. 
4, except for NGC 2244, NGC 2264 and NGC 6611 since for 
them such plots are provided by Sagar (1987). In order to 
convert apparent V values into My we use the relation A„ — 
3.25E{B — V) and distances given in Table 1. We also plot in 
the figure the variation of mean E{B — V) and its standard 
deviation with Mv ■ For estimating these values we grouped 
cluster members in such a way that five or more stars are 
present in a group. An inspection of Fig. 4 indicates that: 

(i) Clusters NGC 663, NGC 869, NGC 884, NGC 1502, 
NGC 1893, Tr 14, Tr 15 and Tr 37 show no correlation of 
colour excess with luminosity. The same behaviour was also 
found by Sagar (1987) for the clusters namely NGC 1976, 
NGC 2244, NGC 2264, NGC 4755 and NGC 6530. 

(ii) In clusters Coll 228 and Tr 16, both mean values and 
scatter in E{B — V) first increase up to Mv ~ -2.0 mag and 
then they decrease with the My of clusters members. The 
rate of increase in E{B — V) with Mv is considerably slower 
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Figure 4. E{B - V) versus Mv diagrams for NGC 663, NGC 
869, NGC 884, NGC 1502, NGC 1893, Tr 14, Coll 228, Tr 15, Tr 
16, Be 86 and Tr 37. The data points are shown in the left panel 
while the right panel shows the plot of mean E{B — V) and its 
standard deviation. 

than the rate of decrease. The correlation between E{B — V) 
and Mv for these clusters is thus complex. 

(iii) In cluster Be 86, there is no variation of E{B — V) 
with Mv for My < 0.0 mag, while for fainter members the 
values of E{B ~ V) decrease with Mv- 

4.4 Variation of colour excess E{B — V) with 
spectral type 

Fig. 5 shows the variation of E{B — V) against spectral class 
of the cluster members. In all the clusters we have only early 
(O, B and A ) type stars. Clusters NGC 663, NGC 869, NGC 
884, NGC 1502, NGC 1893, Tr 14, Tr 15 and Tr 37 show 
no correlation of E{B ~ V) with spectral type. However, for 
Coll 228, Tr 16 and Be 86, E{B — V) decreases as one moves 
from B to A spectral type, while it shows no such variation 
in O type stars. Sagar (1987) also found no corelation of 
E{B — V) with spectral class in the clusters NGC 1976, 
NGC 2244, NGC 2264, NGC 4755, NGC 6530 and NGC 
6611 but in the case of extremely young clusters NGC 6823, 
NGC 6913 and IC 1805, a decreasing trend of E{B - V) 
with spectral class was observed. 

The mean value of E{B — V) and its standard deviation 
with spectral type are also plotted in Fig. 5 for studying the 
variation of scatter in E{B — V) with spectral type. For 
estimating these values we grouped cluster members of one 
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Table 3. Spatial variation of E{B — V) across the cluster Tr 14. The mean values of E(B — V) with their 
standard deviation in mag in 5' x 5' areas are indicated in the appropriate boxes, with the number of stars 
used for this purpose given in brackets. Coordinates are relative to the cluster centre given in Table 2. 



Aa ^ -35 to -30 -30 to -25 -10 to -5 -5 to to 5 5 to 10 

i 



-5 to - - 0.73 0.72±0.3 0.54±0.7 0.47±0.2 

(1) (10) (12) (4) 

to 5 - - 0.82±0.22 0.56±0.9 0.58±0.2 0.44±0.05 

(3) (15) (14) (4) 

10 to 15 0.71±0.10 0.59 - - _ _ 
(4) (1) 



subspectral type. No significant trend is seen between tiie 
scatter in E{B — V) and spectral class. 



4.5 Discussion on non-uniform extinction 

On the basis of Figs. 4 and 5 we can conclude tfiat NGC 
663, NGC 869, NGC 884, NGC 1502, NGC 1893, Tr 14, Tr 
15 and Tr 37 do not show any variation of E{B — V) with 
luminosity as well as spectral class. The clusters Coll 228, 
Tr 16 and Be 86 show a complex correlation of the E{B — V) 
with luminosity as well as spectral class. In these clusters, 
E{B — V) decreses with spectral type and Mv for the B 
stars, while O stars have higher reddening than the B stars. 
The correlation of E{B - V) against My for Tr 16 and Coll 
228 reaches a maximum near Mv = -2.0 mag, then decreases 
for brighter stars. 

The dependency of colour excess E{B — V) with lu- 
minosity and spectral class for the clusters Coll 228, Tr 16 
and Be 86 can be explained with the help of theories of 
star formation. According to Krelowski & Strobel (1981), 
highly reddened luminous stars in association and young 
clusters may have relic circumstellar matter around them. 
On the other hand, theoretical models for star formation 
processes given by Yorke & Krugel (1976) and Bhattachar- 
jee & Williams (1980) suggest that the relative mass of the 
relic envelope present around newly formed massive stars is 
positively correlated with the stellar mass. MS stars emit 
strong stellar winds and UV radiation. The time required 
for them to blow off their relic material will depend upon 
factors like the mass of their relic envelope, resistance of- 
fered by the immediate surroundings, etc. (see Sagar 1987 
for details). Hence one may say that during star formation 
processes in molecular clouds, conditions around MS stars of 
the clusters Coll 228, Tr 16 and Be 86 are such, that relic cir- 
cumstellar material is still present around fainter (M„ > -2.0 
mag) stars, while it has been blown away from the brighter 
stars of these clusters. In addition to this all members in 
clusters where the E{B — V) values are not dependent ei- 
ther on luminosity or spectral class, are either younger or 
have ages comparable to Coll 228, Tr 16 and Be 86. 
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Figure 5. E{B — V) versus spectral class plots for the clusters 
NGC 663, NGC 869, NGC 884, NGC 1502, NGC 1893, Tr 14, 
Coll 228, Tr 15, Tr 16, Be 86 and Tr 37. The data points are 
shown in the left panel while the right panel shows the plot of 
mean E{B — V) and its standard deviation. 



5 EXTINCTION LAW 

In order to understand the properties of interstellar matter 
causing non-uniform extinction in the direction of clusters, 
we plot the ratios of colour excesses E{U ~ V), E{V — J), 
E{V~H) and E{V~K) with E{B-V) against it in Fig. 6 for 
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Figure 6. The colour excess ratios E{U -V)/E(B -V), E{V - J)/E(B -V), E{V - H)/E(B -V) and E{V - K) / E(B -V) are plotted 
against E(B - V) for NGC 663, NGC 869, NGC 884, NGC 1502, NGC 1893, NGC 2244, Tr 14 and Tr 37. The dotted Unes show the 
ratio for a normal interstellar extinction law along with the boundaries expected from observational uncertainties of 30% in the optical 
and 40% in the near-IR (dashes lines). The mean values of the observed colour excess ratios along with their standard deviation and 
number of stars used in deriving them (inside the bracket) are shown on the plot. 



NGC 663, NGC 869, NGC 884, NGC 1502, NGC 1893, NGC 
2244, Tr 14 and Tr 37. The ratio of colour excesses E{U-V), 
E{V-R), E{V-I), E{V-J), E{V-H) and E{V-K) with 
E(B ~ V) for the members of NGC 2264, NGC 6611, Coll 
228, Be 86, Tr 15 and Tr 16 are plotted against E{B - V) 



in Fig. 7. We have considered only those cluster members 
which have E{B ~V) > 0.1 mag, as colour excess ratios for 
stars with smaller E{B — V) values have large uncertainties. 
This criterion has excluded only a few stars in NGC 2264 
since for all other clusters under study E{B — V)min is > 0.1 



© 2001 RAS, MNRAS 000, 



8 R. K. S. Yadav and Ram Sagar 



NGC 3264 NCC 6611 TR 15 Tr 16 Coll 338 Be 86 



Ed 2.0 



21 

- 






"-ll-.. 
























- — 

• 






8.23±6.7a (7) 
. 1 , , , 1 . 


3.23±D.12 (33) 
1 ... 1 . 


|3.Q6±0,47 (ia)| 


. •• 

3.53±0.33 (49) 

1 . . 


- * 3.21± 1.0 (a) 

1 . 1 


2.61±0 15 (10) 

.1 


IS B 


• ■ • 






• • 






• • 


•• • .'• 


= — ■ — i^'Hf- -. 






r *i 








5,28±3-70 (7) 



S,63±0,30 (42) 



2,63±0,57 (14) 

1 


3-a8±0 53 (59) 



3-52±0,33 (6) 

1 


— 5 5 - 


*. '■' ' 


• ••J. 

J.- • • 


• • A 
» •* 


. — 










- ^ - 






3,00±1.35 (7) 

. 1 ... 1 . 


3.02±0 26 (40) 

1 ... 1 


3.57±0,29 (15) 

1 


3.B3±0,30 (51) 

1 . . 


2.40±1,15 (a) 

1 


2.54±0.10 (9) 

.1 


• • • 




• •• 

.. • •• 

V • • . • 


» 


/ 

- • 




1 91±0.90 (7) 

1 ... 1 . 


1.93±0.09 (32) 

1 ... 1 


:.92±D 32 ('l6) 

1 


2.19±0.29 (22) 

1 . . 


2.13±0.14 (9) 

1 , 1 


1.90±0.09 (9) 

.1 


- s : 
1 ' • _ 

^_ • - 


■ . • 




* 

----- r.-- 




- 

— - 
_ 


■ 1.12±0.60 (8) 
, , 1 , , , 1 , 


0.9B±0 06 (33) 
1 ... 1 


O.Q3±0 08 (14) 

1 


0.94±0.11 (22) 

1 , , 1 


Q.99±0.1.4 (10) 

1 , 1 


Q.78±0.05 (7) 

1,1, 


0.2 0.4 






— • • 






1.76±0.17 (10) 

1 ... 1 , 


. • • • 

1.72±0.03 (105) 



. .'••^ 

l.B8±0,29 (16) 

1 


: - 7:- 

- 1,72±0.03 (125) 

, , l' , , 1 


•■ ^•s • • 

1.73±0.03 (55) 

1 


1.72±0,04 (59) - 

1,1, 



, , , \ , LJ , . , I , , , I , I I , , , , , I , . I , LJ I , L_ 

0-2 O-* 0.4 0.8 0.4 0.6 0,6 1.2 0.6 1.2 0.4 0.6 



E(B-V) in mag 



Figure 7. The colour excess ratios E{U-V)/E{B-V}, E{V -R)/ E(B-V), E{V-I)/E(B-V), E{V-J)/E{B-V), E(V -H)/ E(B-V) 
and E{V - K)/E{B - V) are plotted against E(B - V) for NGC 2264, NGC 6611, Tr 15, Tr 16, Coll 228 and Be 86. The dotted and 
short dashed lines and symbols are the same as in Fig. 6. 



mag (see Table 2). In Fig. 6 and 7, horizontal dotted lines 
represent the colour excess ratios corresponding to the nor- 
mal interstellar extinction law described by Mathis (1990). 
Errors expected from the observational uncertainties in the 
colour excess ratios are also shown as short dashed horizon- 
tal line in the figures. They are ~ 30% in the optical but 
become larger, ~ 40%, in the near-IR. The average values 



of the observed colour excess ratios along with its a and 
number of stars used in determining them are also given in 
the appropriate boxes. An inspection of Figs. 6 and 7 clearly 
indicates that:- 

(i) The average value of E{U - V)/E{B - V) for all the 
clusters agree within la, where a is the standard deviation 
of data points, with the normal interstellar value of 1.73. 
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However, there are a few stars in NGC 663, NGC 884, NGC 
1502, NGC 2244, NGC 2264, Tr 14, Tr 15, Tr 16 and Coll 
228 showing ratios that deviate more than their observa- 
tional error, but most are within twice their errors. In Tr 15, 
all colour excess ratios are not too different from the normal 
one except for E{U — V)/E{B — V) where values between 
2.22 and 3.13, compared to the normal value of 1.73, are 
obtained for stars with E{B-V) > 0.45. Tapia et al. (1988) 
also noticed this anomaly. The least square linear regression 
to the data points yields '^[b-v] ^ (3.82±0.88)-B(B - V) 
+ (0.27±0.33). Tapia et al. (1988) conclude that Tr 15 is lo- 
cated in a denser intracluster dust cloud whose grains have 
been 'processed' in such a way as to produce disproportion- 
ally larger E{U — V) excess per unit column density as com- 
pared to the average interstellar medium. In all the clusters, 
attenuating material seems to have the properties of normal 
interstellar matter except for a few stars whose anomalous 
extinction arises from their being a field star or a cluster 
star with peculiar intrinsic colours. 

(ii) In the clusters NGC 663, NGC 869, NGC 884 and 
NGC 1502 the colour excess ratios for members at A > Aj 
are generally smaller than the normal values, when com- 
pared with E{B — V) (see Fig. 6). This and the values of R 
{=1.1E{V - K)/E{B - V) following Whittet & van Breda 
(1980)) of these clusters indicate that the size distribution of 
the dust particles in the direction of these clusters is biased 
towards smaller values than the size of normal interstellar 
matter such that they produce smaller values of colour ex- 
cess ratios when compared with E{B — V), but normal values 
while compared with E{V — J). Creze (1972) also indicates 
that in the direction of all these clusters, the value of R is 
relatively lower than other sight lines through the Galaxy. 
Tapia et al. (1991) have also shown that the size distribution 
of dust particles in the direction of the cluster NGC 1502 
is biased towards smaller values. All these clusters have a 
similar age (~ 12 Myr) and are located in a small sector of 
the Galactic plane with 130° < 1 < 145°. We therefore con- 
clude that the interstellar matter in this part of the Galaxy 
is somewhat different from the normal one. 

(iii) The colour excess ratios of the members in NGC 
1893, NGC 2244, Tr 37 and Be 86 seems to be generally 
normal at all wavelengths from U to K passbands (see Fig. 
6 and 7). This indicates that the extinction law is close to 
normal galactic extinction in the direction of these clusters. 
Perez et al. (1987) and MorbideUi et al. (1997 ) have also 
shown that the value of R is normal in the direction of NGC 
2244 and Tr 37 respectively. In the direction of NGC 1893, 
Tapia et al. (1991) have also indicated an average galactic 
extinction law. 

(iv) The colour excess ratios to E{B — V) of the cluster 
members are larger than the normal values at A > Xr in 
NGC 6611, Coll 228 and Tr 16; and at A > Aj in Tr 14, 
but the values become normal when compared to E{V — J) 
except in NGC 6611 where it is slightly lower than the corre- 
sponding normal values. This implies that the extinction law 
is anomalous in the sense that, in these regions, grain size 
distribution seems to be biased towards larger than normal 
size particles. Hillenbrand et al. (1993) have also reported a 
higher value of R for NGC 6611. 

(v) The cluster NGC 2264 exhibits a large scatter from 
the normal values of the colour excess ratios. Perez et 
al. (1987) also found a large deviation in the value of R 
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Figure 8. Plots of near-IR flux excess/deficiency in terms of A(V- 
H) and A(V-K) against the colour excess E{V — J) (see text). The 
horizontal dotted lines denote zero excess. The short dashed lines 
denote the extent of the expected errors. 



most probably due to the presence of circumstellar material 
around some cluster members (see next section). 



6 NEAR-IR FLUXES 

Contributions towards the reddening measured in terms of 
colour excesses of the cluster members include (i) interstellar 
extinction between the observer and the cluster, (ii) intra- 
cluster extinction due to the gas and dust present inside 
the cluster; and (iii) circumstellar extinction due to rem- 
nant disks or envelopes, which have not got fully dissipated 
from the cluster members because of their youth. If one as- 
sumes that the interstellar material in the foreground of the 
clusters follows the normal interstellar extinction law, then 
whatever deviations are seen from the normal values in the 
colour excess ratios of stars is only due to (ii) or (iii) or 
both. The colour excess used to measure interstellar matter 
should, however be selected carefully as the clusters under 
study contain young stellar objects. In such circumstances 
one uses E{V-J) instead of E{B-V) (Smith 1987; Tapia et 
al. 1988,1991; Sagar and Qian 1989, 1990), because E{V-J) 
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Table 4. Near-IR flux excess/deficiency in members of clusters under study. 



cluster Webda Star No. V E{V - J) A(V-H) A(V-K) Sp Type 

(mag) (mag) (mag) (mag) 



NGC 884 


1781 


9.26 


1.57 


-0.60 


-0.58 


BlIV 


NGC 2264 


90 


12.71 


1.17 


0.91 


1.98 


B4V 




100 


10.04 


0.64 


0.93 


2.16 


A2IV 




165 


10.99 


0.49 


0.43 


1.00 


A2V 




46 


9.18 


0.70 


0.09 


0.72 


A3V 


Tr 14 


15 


11.85 


1.88 


0.53 


1.27 


B7V 


Tr 16 


68 


12.14 


1.50 


0.63 


1.71 


05V 



does not depend on environmental properties such as chem- 
ical composition, shape, structure and the degree of align- 
ment of the interstellar matter (Vashchinnikov & Il'in 1987; 
Cardelli, Clayton & Mathis 1989). The differences between 
the observed colour excess in (V— H) and (V— K) based on 
spectral classification and the derived colour excess from 
E{V — J) assuming normal extinction law are calculated. 
These differences are plotted with E{V — J) in Fig. 8. The 
major sources of error in the differences are the observational 
uncertainties in JffK magnitudes, inaccuracies in estimation 
of E(y — J) and errors in the spectral classification. The dif- 
ferences can therefore be considered statistically significant 
only if their absolute values are larger than ~0.5 mag. 

An inspection of Fig. 8 indicates that the values of 
A(V— H) and A(V— K) are distributed randomly around the 
horizontal line, except in the case of NGC 6611, where they 
are systematically lower. Most of the cluster members have 
absolute values of A(V— H) and A(V— K) close to zero in- 
dicating that there are no near-IR excess fluxes. However, 
a few stars show an excess and deficiency in radiation at 
H and K wavelengths. They are listed in Table 4. We can 
clearly see from this table that all the stars are producing 
excess fiux more in K compared to H, except one star NGC 
884-1781 which has approximately equal fiux deficiency in 
both bands. Strom et al. (1971, 1972) noticed circumstel- 
lar envelope around star NGC 2264-90. Sung, Bessell & Lee 
(1997) identified NGC 2264-100 as a pre- main sequence star, 
while Young (1978) found emission in star NGC 2264- 
64. Smith (1987) found an infrared excess in star Tr 14-15. 
These observations, thus support our results listed in Table 
4. We have identified for the first time the presence of a 
near- infrared excesses in NGC 2264-165 and Tr 16-68 and a 
deficiency in NGC 884-1781. The stars are of O, B and A 
type main sequence, except NGC 884-1781 and NGC 2264- 
100. 



7 CONCLUSIONS 

From the present analysis of reddening and ratios of colour 
excesses of members in 14 young (age < 20 Myr) open clus- 
ters, the following main conclusions can be drawn: 

(i) Non-uniform extinction decreases systematically with 
the age of the clusters. We therefore conclude that young 
open star clusters (age < 10^ yr) still have unused gas and 
dust inside them, while in older clusters (age > 10* yr) , they 
are either used in star formation or blown away by the UV 
radiation and strong winds of hot stars or both. 



(u) The cluster members of Coll 228, Tr 16 and Be 86 
show significant variation of E{B — V) with luminosity and 
spectral class. This suggests that some stars have retained 
their circumstellar material whilst it has been stripped away 
from most others. In these clusters, the colour excess ratios 
at A > Xr have larger values than normal. Ifowever, the 
presence of near-IR excess is not observed at A < 2.2 /im in 
most of the members of Coll 228, Tr 16 and Be 86 except for 
one star (see Table 4). Observations at longer wavelengths 
are required to confirm the presence or absence of circum- 
stellar dust. 

(iii) There is no significant variation of colour excess spa- 
tially except in Tr 14. In the clusters NGC 663, 869, 884, 
1893, Tr 16 and Be 86, random variation of E{B - V) is 
present. In the southern part of the cluster Tr 14, the value 
of E{B — V) increases systematically as one moves from east 
to west (see Table 3). 

(iv) Colour excess diagrams show that cluster members 
generally follow the normal interstellar extinction law at 
optical wavelengths, while anomalous behaviour has been 
noticed at A > Aj for stars in some clusters (see Figs 6 and 
7). In the direction of a small part of the Galaxy ( longitude 
ranging from ~ 130 to 145 deg along the Galactic plane) 
the size distribution of the dust particles is biased towards 
smaller values than normal. 

(v) Some stars show anomalous behaviour at H and K. 
The near-IR excess radiation, which generally increases with 
wavelengths from H to K could be due to the presence of hot 
dust close to the stars or perhaps a companion? In contrast, 
star NGC 884-1781 show the same near-IR flux deflciency 
in both H and K passband. 

On the basis of the above, one may conclude that there is no 
uniformity in extinction properties amongst these clusters. 
In fact, they differ from cluster to cluster. All this indicates 
that non-uniform extinction observed in young clusters can- 
not be understood in terms of a simple physical scenario. 
Actually, it depends upon factors such as the age of the 
cluster members, initial spatial distribution of matter in the 
molecular clouds, sequential star formation processes and 
the distribution of hot O and B stars in the cluster etc. A 
complicated physical scenario is therefore required to ex- 
plain the non-uniform extinction in young star clusters. 
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